Treatment with 2 commercial antibiotics reduced clinical and systemic signs of pneumonia and the abundance of pathogenic bacteria in the upper respiratory tract of preweaning dairy calves.

The aim of this study was to evaluate the effect of therapeutically administered tildipirosin or florfenicol + flunixin meglumine for the treatment of bovine respiratory disease (BRD) accompanied by fever in calves before weaning compared with diseased and untreated animals. As specific objectives, we evaluated the composition of the bacterial microbiota of the upper respiratory tract (URT) and blood and health parameters of the animals. Preweaning Holstein female calves diagnosed with naturally acquired pneumonia were randomly assigned to one of the following experimental groups on the day of diagnosis (d 0): (1) TLD (n = 36): single subcutaneous injection with 4 mg/kg tildipirosin; (2) FLF (n = 33): single subcutaneous injection with an antimicrobial (40 mg/kg florfenicol) combined with a nonsteroidal anti-inflammatory drug (2.2 mg/kg flunixin meglumine); and (3) NEG (n = 35): no treatment within the first 5 d following enrollment. The NEG treatment group was closely monitored for 5 d, and calves were removed from the study following a standardized late treatment protocol, when necessary, to minimize health concerns. Healthy untreated calves (CTR; n = 31) were also selected for the study and used as controls. Blood samples used for biochemical analysis and nasopharyngeal swabs used for evaluation of URT microbiota were collected daily from d 0 until d 5 and then weekly until weaning. Next-generation sequencing of the 16S rRNA gene was used to assess the URT microbiota at the phylum and genus levels. Clinical signs associated with pneumonia and otitis media were assessed daily, as was the need for antibiotic interventions. Calves in the TLD and FLF groups had faster recovery from fever within the first 5 d after enrollment. In addition, antibiotic-treated calves reached the same serum haptoglobin levels as healthy calves on d 2 after diagnosis, whereas calves in the NEG group had higher haptoglobin levels than the CTR group until at least d 5 after BRD diagnosis. Calves in the TLD and FLF groups had a lower risk of treatment for pneumonia (FLF = 22.8%; TLD = 27.7%) from d 5 to weaning than calves in the NEG group (54.7%). Furthermore, FLF treatment had a significantly lower risk of nasal discharge, otitis media, and treatment failure compared with the NEG group, but did not differ from the TLD group. Differences in the composition of the URT microbiota were found between groups, and the genus Mycoplasma was the most abundant in samples collected from the URT of calves with and without pneumonia. Both drugs were effective in reducing the mean relative abundance (MRA) of important genera associated with pneumonia (Mannheimia and Pasteurella), although an increase in Mycoplasma MRA was observed for tildipirosin-treated calves. In conclusion, both drugs were effective in reducing the inflammatory signs of pneumonia and the need for antimicrobial treatment after enrollment compared with no treatment. In addition, both TLD and FLF were effective in reducing the MRA of important bacterial genera associated with pneumonia; however, TLD treatment was associated with increased Mycoplasma MRA compared with healthy and untreated calves.

Heat stress response in slick vs normal-haired Criollo Limonero heifers in a tropical environment.

To assess the effect of hair type on the heat stress response, 20 Criollo Limonero heifers with slick (n = 11) or normal hair (n = 9) were studied. Under a high temperature-humidity index (THI) environment, heat stress response was assessed through physiological variables that included respiration rate (RR), heart rate (HR), ruminal frequency (RMF), rectal temperature (RT), saliva pH (SPH), and lymphocyte count (LC) in the morning (5:00 AM, 27.4 °C, 64% relative humidity, THI = 77) and afternoon (1:00 PM, 34.5 °C, 70% relative humidity, THI = 88). A case-control study using a split plot design was used. Data were analyzed using ANOVA (PROC MIXED SAS 2010) and a statistical model comprising the fixed effects of hair length, sampling hour, interaction of hair length by sampling hour, and the random effect of animal nested within hair type on physiological variables associated with heat stress response. Sampling hour influenced (P < 0.0001) RR, RT, and (P < 0.003) SPH. Hair length influenced RR (P < 0.01) and RT (P < 0.04) and tended to influence LC (P < 0.07). The interaction of sampling hour by hair influenced RR (P < 0.04), RT (P < 0.0002), and both SPH and LC (P < 0.05). During afternoon hours, slick-haired heifers had lower values for RR (81 ± 4.2 vs 102 ± 4.7 bpm; P < 0.01), RT (39.5 ± 0.1 vs 40.3 ± 0.1 C°; P < 0.002), and LC (60 ± 3.2 vs 72.3 ± 3.6; P < 0.09) than normal-haired heifers. In normal-haired heifers, SPH increased during afternoon compared to morning-hours (8.66 ± 0.1 vs 9.11 ± 0.1; P < 0.04). It was concluded that slick-coated heifers exhibited an enhanced capability to cope with heat stress compared to normal-haired heifers likely due to an enhanced capacity for heat dissipation.

All sheeps and sizes: a genetic investigation of mature body size across sheep breeds reveals a polygenic nature.

Mature body size is genetically correlated with growth rate, an important economic trait in the sheep industry. Mature body size has been studied extensively in humans as well as cattle and other domestic animal populations but not in sheep. Six-hundred and sixteen ewes, across 22 breeds, were measured for 28 linear measurements representing various skeletal parts. PCA from these measures generated principal components 1 and 2 which represented 66 and 7% of the phenotypic variation respectively. Two-hundred and twenty sheep were genotyped on the Illumina Ovine HD beadchip for a GWAS investigating mature body size and linear body measurements. Forty-six (Bonferroni P < 0.05) SNP associations across 14 chromosomes were identified utilizing principal component 1, representing overall body size, revealing mature body size to have fewer loci of large effect than other domestic species such as dogs and horses. Genome-wide associations for individual linear measures identified major quantitative trait loci for withers height and ear length. Withers height was associated (Bonferroni P < 0.05) with 12 SNPs across six chromosomes whereas ear length was associated with a single locus on chromosome 3, containing MSRB3. This analysis identified several loci known to be associated with mature body size in other species such as NCAPG, LCORL, and HMGA2. Mature body size is more polygenic in sheep than other domesticated species, making the development of genomic selection for the trait the most efficient option for maintaining or reducing mature body size in sheep.

Exploring physiological and genetic variation of digital cushion thickness in Holstein and Jersey cows and bulls.

The objectives of our study were to expand phenotypic characterization of digital cushion thickness (DCT) to the Jersey breed and include mature bulls and to identify breed-specific quantitative trait loci (QTL) for DCT within the Jersey or Holstein breeds and common QTL across breeds and sexes to better understand biological regulation and guide future use of marker-assisted selection. In a cohort of 698 cows and 85 bulls (Holstein and Jersey) from 8 farms in New York State, data were collected on DCT, body condition score (BCS), withers height, and sacral height. All animals underwent sonographic examination of the digital cushion evaluated at the sole ulcer site for the right front and hind feet. Linear mixed models were conducted on DCT separately for cows and bulls with fixed effects of time point, breed, age, digit, and BCS group. The models included random effects to control for the random subset of animals per farm, repeated measures, and multiple measurements from each animal. The phenotypic results indicated that DCT varied by time point, breed, age, digit, and BCS group for cows and by breed, age, digit, and BCS group for bulls. For the genotypic study, 616 cow DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA), whereas 76 bull DNA samples were genotyped on different platforms ranging from 5K to 150K. Multiple genome-wide association studies were conducted to highlight pertinent phenotyping parameters and genetic markers for genomic selection. Data were separated into 8 data sets based on different combinations of breed and sex. Each data set was assessed for quality of markers and samples before conducting genome-wide association studies for DCT, testing the inheritance models and genetic variation of digit, foot, and average thickness. Ten markers passed the Bonferroni correction threshold and 9 passed false discovery rate from 10 genome-wide association studies using a combination of the covariates breed, sex, genotyping batch plate, age, BCS, withers height, and sacral height. Of the 43 candidate genes, 8 novel biologically plausible genes were identified on Bos taurus autosomes 3, 4, 7, and 9: SFRS18 and LRRFIP1 function in fat deposition, whereas AHR, BZW2, EFNA5, USP45, and VAV3 effect bone growth, and SOSTDC1 is related to epidermal keratinocyte function. The genetic markers associated with DCT in this study were explored for variation between cows and bulls within and across breeds for their potential use in marker-assisted selection.

Phenotypic characterization and genome-wide association studies of digital cushion thickness in Holstein cows.

The bovine digital cushion is a compression pad between the distal phalanx and sole and has been associated with claw horn disruption lesions. Digital cushion thickness (DCT) is estimated to be moderately heritable. Therefore, the objectives of our study were to examine influences of management and environment on DCT and to identify genetic markers and candidate genes associated with DCT. In a cohort of 502 Holsteins from 5 farms in New York State, DCT and body condition score (BCS) were collected twice, at <137 d prepartum and from 86 to 127 d in milk, corresponding to periods when the digital cushion is thickest and thinnest, respectively, as determined by previous research. Cows underwent sonographic examination of the digital cushion evaluated at the typical sole ulcer site for the right front and hind foot. Linear mixed models were conducted on DCT with the fixed effects of time point, digit, wither height, sacral height, BCS group, and multiple farm system variables separately and included random effects to control for the random subset of cows per farm, repeated measures, and multiple measurements from each cow. The phenotypic results indicated that DCT varied by sample time point, sacral height, parity, digit, BCS group, and wither height. For the genotypic study, 447 DNA samples were genotyped on the Illumina BovineHD 777K BeadChip (Illumina Inc., San Diego, CA). Quality assessment of markers and samples provided a final data set of 431 samples and 579,449 markers. Genome-wide association studies were conducted for DCT testing inheritance models and genetic variation of digit, foot, time point, and average thickness. One marker passed the Bonferroni correction threshold and 26 passed false discovery rate from 4 genome-wide association studies with covariates of sequencing batch plate, parity group, BCS, wither height, and sacral height. Ten candidate genes were identified, with 2 genes on Bos taurus autosomes 24 and 29 involved in biological functions related to the digital cushion: MC4R and DLG2 were related to fat deposition and bone growth, respectively. The genetic markers discovered in this study have the opportunity to be used in breeding programs using genomic selection to select against claw horn disruption lesions and lameness due to associations between the markers and DCT. Further studies on the biologically plausible candidate genes may identify causative genetic variants and how they relate to DCT through gene regulation, expression, structure, or copy number variation.

Population genomic structure and linkage disequilibrium analysis of South African goat breeds using genome-wide SNP data.

The sustainability of goat farming in marginal areas of southern Africa depends on local breeds that are adapted to specific agro-ecological conditions. Unimproved non-descript goats are the main genetic resources used for the development of commercial meat-type breeds of South Africa. Little is known about genetic diversity and the genetics of adaptation of these indigenous goat populations. This study investigated the genetic diversity, population structure and breed relations, linkage disequilibrium, effective population size and persistence of gametic phase in goat populations of South Africa. Three locally developed meat-type breeds of the Boer (n = 33), Savanna (n = 31), Kalahari Red (n = 40), a feral breed of Tankwa (n = 25) and unimproved non-descript village ecotypes (n = 110) from four goat-producing provinces of the Eastern Cape, KwaZulu-Natal, Limpopo and North West were assessed using the Illumina Goat 50K SNP Bead Chip assay. The proportion of SNPs with minor allele frequencies >0.05 ranged from 84.22% in the Tankwa to 97.58% in the Xhosa ecotype, with a mean of 0.32 ± 0.13 across populations. Principal components analysis, admixture and pairwise FST identified Tankwa as a genetically distinct population and supported clustering of the populations according to their historical origins. Genome-wide FST identified 101 markers potentially under positive selection in the Tankwa. Average linkage disequilibrium was highest in the Tankwa (r(2)  = 0.25 ± 0.26) and lowest in the village ecotypes (r(2) range = 0.09 ± 0.12 to 0.11 ± 0.14). We observed an effective population size of <150 for all populations 13 generations ago. The estimated correlations for all breed pairs were lower than 0.80 at marker distances >100 kb with the exception of those in Savanna and Tswana populations. This study highlights the high level of genetic diversity in South African indigenous goats as well as the utility of the genome-wide SNP marker panels in genetic studies of these populations.

Genomic evaluation, breed identification, and population structure of Guernsey cattle in North America, Great Britain, and the Isle of Guernsey.

As of December 2015, 2,376 Guernsey bulls and cows had genotypes from collaboration between the United States, Canada, the United Kingdom, and the Isle of Guernsey. Of those, 439 bulls and 504 cows had traditional US evaluations, which provided sufficient data to justify investigation of the possible benefits of genomic evaluation for the Guernsey breed. Evaluation accuracy was assessed using a traditional 4-yr cutoff study. Twenty-two traits were analyzed (5 yield traits, 3 functional traits, and 14 conformation traits). Mean reliability gain over that for parent average was 16.8 percentage points across traits, which compares with 8.2, 18.5, 20.0, and 32.6 percentage points reported for Ayrshires, Brown Swiss, Jerseys, and Holsteins, respectively. Highest Guernsey reliability gains were for rump width (44.5 percentage points) and dairy form (40.5 percentage points); lowest gains were for teat length (1.9 percentage points) and rear legs (side view) (2.3 percentage points). Slight reliability losses (1.5 to 4.5 percentage points) were found for udder cleft, final score, and udder depth as well as a larger loss (13.6 percentage points) for fore udder attachment. Twenty-one single nucleotide polymorphisms were identified for Guernsey breed determination and can be used in routine genotype quality control to confirm breed and identify crossbreds. No haplotypes that affect fertility were identified from the current data set. Principal component analysis showed some divergence of US and Isle of Guernsey subpopulations. However, the overlap of US, Canadian, UK, and Isle of Guernsey subpopulations indicated the presence of gene flow, and the similarities in the subpopulations supports a common genomic evaluation system across the regions.

The SLICK hair locus derived from Senepol cattle confers thermotolerance to intensively managed lactating Holstein cows.

The SLICK haplotype (http://omia.angis.org.au/OMIA001372/9913/) in cattle confers animals with a short and sleek hair coat. Originally identified in Senepol cattle, the gene has been introduced into Holsteins. The objectives of the current study were to determine (1) whether lactating Holsteins with the slick hair phenotype have superior ability for thermoregulation compared with wild-type cows or relatives not inheriting the SLICK haplotype, and (2) whether seasonal depression in milk yield would be reduced in SLICK cows. In experiment 1, diurnal variation in vaginal temperature in the summer was monitored for cows housed in a freestall barn with fans and sprinklers. Vaginal temperatures were lower in slick-haired cows than in relatives and wild-type cows. In experiment 2, acute responses to heat stress were monitored after cows were moved to a dry lot in which the only heat abatement was shade cloth. The increases in rectal temperature and respiration rate caused by heat stress during the day were lower for slick cows than for relatives or wild-type cows. Moreover, sweating rate was higher for slick cows than for cows of the other 2 types. In experiment 3, effects of season of calving (summer vs. winter) on milk yield and composition were determined. Compared with milk yield of cows calving in winter, milk yield during the first 90 d in milk was lower for cows calving in the summer. However, this reduction was less pronounced for slick cows than for wild-type cows. In conclusion, Holsteins with slick hair have superior thermoregulatory ability compared with non-slick animals and experience a less drastic depression in milk yield during the summer.

Association of hypothyroid disease in Doberman Pinscher dogs with a rare major histocompatibility complex DLA class II haplotype.

Canine hypothyroid disease is similar to Hashimoto's disease in humans, which has been shown to be associated with human major histocompatibility complex (MHC) genes. We have collected 27 Doberman Pinschers affected with primary hypothyroid disease and compared their MHC class II haplotypes with 129 unaffected Doberman Pinschers. Three dog-leucocyte antigen (DLA) genes, DLA-DRB1, DQA1 and DQB1, were characterized by sequence-based typing and assigned to haplotypes for each dog. One rare haplotype was found at an increased frequency in the affected dogs compared to the unaffected dogs (Odds ratio = 2.43, P < 0.02). This haplotype has only been found in Doberman Pinschers and Labradors to date.